A circuit QED architecture with current-biased flux qubits 
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We theoretically study a circuit quantum electrodynamics (QED) architecture with current-biased 
flux qubits. The qubit is coupled to the transmission line resonator by a bias current originating 
from the current mode of the resonator. Ultrastrong coupling regime can be obtained by varying the 
capacitance between the qubit and the resonator. We propose a scalable design for the circuit QED 
with current-biased flux qubits, where the dc-SQUIDs take the role of switching the qubit-resonator 
coupling. An exact calculation on two-qubit coupling strength in the scalable design shows the 
transition from ferromagnetic to antiferromagnetic xy-type interaction. 

PACS numbers: 74.50.-|-r, 03.67.Lx, 85.25.Cp 



I. INTRODUCTION 

An artificial two level system can be coupled with 
the quantized electromagnetic field in a superconduct- 
ing transmission line resonator, while the natural atom 
is coupled with 3D cavity. This circuit quantum electro- 
dynamics (QED) architecturei*^ is a solid-state analog of 
3D cavity QED, providing a strong coupling strength be- 
tween the qubit and resonator owing to the large dipole 
moment of the artificial qubit. The circuit QED scheme 
has been applied to the superconducting qubits. Among 
the superconducting qubits the fiux qubilj^"— has the ad- 
vantage of fast gate operation because the flux qubit does 
not require low anharmonicity in circuit QED scheme. 
There have been many studies for the circuit QED with 
the superconducting flux qubit.— However, the induc- 
tive coupling between the flux qubit and the transmission 
line resonator of the circuit is too weak to perform the 
quantum gate operation. 

Recently a galvanic coupling scheme for the circuit 
QED with the flux qubits has been proposed to en- 
hance the coupling strength by sharing the flux qubit 
loop and the resonator transmission line.—"— In this study, 
we propose a bias current coupling scheme between the 
flux qubit and the transmission line resonator; the qubit 
and the resonator are not galvanically coupled with each 
other, but by a current flowing through the capacitance 
between the qubit and the resonator. In this scheme 
the three-junctions flux qubit is coupled by a bias cur- 
rent similarly to the superconducting phase qubit.—"— 
While the states of phase qubit are defined in terms of the 
phase degrees of freedom in a washboard type potential, 
the present current-biased flux qubit uses the persistent 
current states as qubit states. The qubit state prepa- 
ration and the quantum gate operation are achieved by 
the bias current. Our qubit thus has the advantages of 
phase qubit such as fast qubit operation and readout and 
individual addressing. 

In the design of our qubit the capacitance can be con- 
trolled by varying the width of capacitance line extended 
from the qubit loop and the distance between the capac- 
itance line and the transmission line resonator. The cou- 



pling strength between the qubit and the resonator shows 
a maximum where the coupling constant can reach even 
ultrastrong coupling regime with reasonable parameter 
values. 

We also introduce a scalable design with switching 
function and study the behavior of coupling strength for 
various number of qubits. The coupling between qubit 
and the resonator can be switched on/off by using a dc- 
SQUID inserted between the flux qubit and the ground 
plane of circuit QED. Further we analyze the xy-type in- 
teraction between two qubits. The two-qubit coupling 
is shown also to be strong, which requires an exact, not 
perturbative, representation of the two-qubit interaction. 
The obtained two-qubit coupling shows that we can se- 
lectively choose between ferromagnetic and antiferromag- 
netic xy-type coupling by varying system parameters. 



II. A STRONG COUPLING AND SINGLE 
QUBIT OPERATION 

A. Strong qubit-resonator coupling 

Usually the transmon qubit is coupled with the volt- 
age mode of the transmission line resonator through a 
capacitance^^ For superconducting flux qubits, there 
have been many studies to couple the flux qubit with 
the transmission line resonator by using mutual induc- 
tance between the qubit loop and the resonatoriSiiS or 
by sharing the qubit loop with the resonator<^— On the 
other hand the three-junctions flux qubit can also be cou- 
pled with the transmission line resonator through a ca- 
pacitance, but in this case it is coupled with the current 
mode of the resonator. An oscillating bias current flowing 
across the capacitance gives rise to the coupling between 
the qubit and the resonator 

If the three-junctions flux qubit is penetrated by a 
magnetic flux of half flux quantum $0/2, there are two 
current states in the qubit loop. The clockwise and coun- 
terclockwise current states correspond to the local min- 
ima in the effective potential of the qubit loop. If the 
three Josephson junctions are identical, they have equal 
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phase difference ±a with a = 7r/3, depending on the cur- 
rent states. When a bias current /q is apphed to the flux 
qubit where three Josephson junctions are located asym- 
metrically in the loop as shown in Fig. [TJ the coupling 
strength is given by a product of flux variable and bias 
curreniji^ 



ZTT 



(1) 



similarly to the superconducting phase qubit J^i^^ Here 
we have the phase difference a instead of So; because 
the two phases in different sides of the qubit are can- 
celled out each other. This bias current coupling has re- 
cently been realized in experiments with three-junctions 
flux qubit 

In this study, we consider a qubit design which can 
control the capacitance between the qubit and the res- 
onator. Large capacitance allows high bias current, and 
thus makes the qubit-resonator coupling reach even ul- 
trastrong coupling regime. We consider a qubit design 
shown in Fig. [ija), where the three-junctions flux qubit 
is coupled with the transmission line resonator by a ca- 
pacitance line extended from the qubit loop. The width 
w of the capacitance line and the distance d between the 
capacitance line and the resonator can be adjusted to de- 
termine the capacitance between the qubit and the res- 
onator. In this case the capacitance density around the 
qubit location is larger than other area of the resonator, 
providing a strong qubit-resonator coupling. Two capac- 
itors at the ends of the resonator are introduced for the 



current mode of the resonator to be periodic in a scalable 
design. Here we use the second harmonics of the current 
mode, and the arrows show the corresponding current in 
the circuit. 

The microwave passing through the uniform resonator 
in the circuit QED architecture can be described as a 
one-dimensional motion by the Lagrangian 



c{e,e-t) 



l{^6{x,t)f-^^{^e{x,t)f]dx, (2) 



where / and c are the inductance and the capacitance per 
unit length of the uniform transmission line resonator, re- 
spectively, and 9{x,t) = /^^/2 dx'q{x',t) with the linear 
charge density q{x) is a collective field variable. 

Figure [TJb) shows the schematic diagram of (a), where 
c and c' are the capacitance density between the res- 
onator and the ground plane and between the resonator 
and the qubit, respectively. When c' ^ c, almost current 
flows through the qubit at the center and the capacitors 
at the ends of the resonator. The resonator thus is not 
uniform any more, and the equation of motion of the 
field variable in the sector i of the resonator is given by 
in terms of the Euler-Lagrange equation 



1 d^e,{x,t) 



(3) 



If the field variable is represented as a product of spa- 
tial and temporal parts 9i{x,t) = Xi(x)(l)(t), we get the 
equation 
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FIG. 1: (a) A current-biased three-junctions flux qubit where 
Wo is the width of qubit loop and do the distance between the 
resonator and the ground plane. The upper superconducting 
plane is the transmission line resonator and the lower is the 
ground plane, (b) A schematic diagram for the circuit in (a), 
where c and c are the capacitance density. The three Joseph- 
son junctions with phase difference are located asymmet- 
ricaUy along the qubit loop with threading flux $. 
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Ci dx^ 



X^{x) + iLolX^ix) = 0. 



(4) 



The spatial part Xi{x) (—2 < z < 2) is explicitly written 
as 



X{x)^< 



I ^26'^^ 



+ S_2e-*^^ (-f < x < -f + f 
S_ie"*T^^ (-■§ -I- f < a; ^ - 2 



-Boe-*¥^ (-f<x<f) 
-Sie-*^^ (f<x<f-f) 



(5) 



and (jjr ~ —h= ^ is the resonator frequency. Here Ci — c 



and ji — ji for odd i and c; = c' and ji — j2 for even i, 
and thus we have 



1 jlTT 



1 j27r 



(6) 



From Eq. (j4|) we can readily observe that 
{l/ci)dXi{x)/dx and Xi{x) are continuous at the bound- 
ary between the sectors, which means the continu- 
ity of electric potential Vi{x,t) = {l/ci)dO{x,t)/dx = 
(l/c,)VX,;(a;) and current I{x,t) = de{x,t)/dt = 
X{x)(j){t) at boundary. In Fig. [T]the resonator and the 
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qubit are coupled by a bias current flowing into the qubit 
through the capacitance in the region — w/2 < x < wjl. 
The bias current is given by 



w/2 



w/2 



q(x, t)dx 



'If'' 



(7) 



j2 are integers, but in general they are rational numbers 
depending on the ratio c'/c. 

The set (i) of boundary conditions in Appendix A pro- 
vides the relation A-i = —Bi, and the coefficients Ai and 
Bi are determined from the set (ii) of boundary condi- 
tions as 



from the current conservation condition q{x, t) = 
dl{x,t)/dx in the resonator. The other current flows 
from the resonator to the ground plane directly through 
the capacitors with small capacitance density c in Fig. 

HTb). 

In Appendix A we present the boundary conditions 
for general N qubit case. We, first of all, consider N = 1 
case. The boundary conditions are categorized into two 
sets in Appendix A. X{x) should be an odd function 
in order to have a finite coupling constant g, because 
the bias current of Eq. (O vanishes if X{x) is an even 
function for —wjl < x < w/2. In Appendix A, thus, 
the determinant of matrix corresponding to set (ii) of 
boundary conditions should be zero, resulting in 



Cj2 cos — Y 



ic'ji sin^f 



CJ2 COS^f 



ic'jism^-^ 



(8) 



The values of ji, j2 and uJr arc determined from Eqs. ^ 
and ([8]). For uniform capacitance density c! — c, ji and 
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FIG. 2: (a) Current profile of single qubit case for various 
w/wo with do/d = 10 and wq/L = 10"''. The gap S around 
x/L — increases along with w/wq, which gives rise to a 
strong coupling constant g. (b) The enlarged current profile 
around the gap inside the dotted ellipse in (a) for various 
w/wo with do/d = 10 and (c) for various do/d with u;/wo=20. 
Thin dotted lines indicate the boundary position x = w/2, the 
end of the capacitance line. As w/wq and do/d increase, the 
gap grows and is finally saturated. 
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The Lagrangian in Eq. ([2]) of the resonator modes 
can be written as £(0, 4>) = L — ■^^4'^^ with di- 

mensionless constant ^ = {^/L)J2iJ_i^/2-^ii^)'^^ ^'^'^ 
K = {l/L)J2^J^L%(^/c^){VXi{x))^dx. If we introduce 
the representations — 



2^ 

1 / hUrC 



-(a - a'). 



/2k 



(a-fa'f). 



(10) 



(11) 



the Hamiltonian of the resonator modes is written in 
a diagonalizcd form = fujjr{a^a + ^)- The current 
I{x,t) = X{x)(j){t) is then given by 



/(x,t) 



Xix) 



(12) 



where the remaining coefficient A2 in Eq. ([9]) is a com- 
mon factor in the numerator and denominator, and thus 
is cancelled out. 

The total Hamiltonian Hjc ~ Hr + Hq + Hj given by 
the sum of the Hamiltonian for the resonator modes, for 
the qubit, and for the interaction between the resonator 
modes and the qubit is written as a Jaynes-Cummings 
type Hamiltoniani^ 

Hjc ^ fi^^ra^a+'^o-z + igc^xia- a^), (13) 

where Wa is the qubit frequency and the last term rep- 
resents the coupling between the qubit and the current 
mode in the resonator, which is different from gax{a + a^) 
in transmon case. This type of bias current coupling term 
for the three-junctions flux qubit has also been derived 
recently in a different manner 

The amplitude of bias current in Eq. (O is given by 



In = 



U 



(14) 
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with 



(5 = 2 



X(f) 



(15) 



and thus the couphng strength g in Eq. ([T|) is 
determined by 5 and Ur- Since the amplitude of 
current I{x,t) in Eq. (fT2|) satisfies the condition, 

i^/L)J2^I-i%iX^ix)/V^)^dx = 0.5, ^ has the maxi- 
mum value of \/2 when the current profile takes the rect- 
angular function form. 

In Fig. [5] (a) we show the current profile I{x) = 

where we set Wq/L = 10"'' corresponding 
to Wo = l^m when L = 10mm. Here a finite current gap 
develops around the qubit location at the center of the 
resonator. The capacitance c' = (wdo/wod)c between 
the resonator and the capacitance line increases along 
with w/wo and do/d, which enables more charges to flow 
across the qubit. The resulting large bias current, pro- 
ducing the gap in current profile of Fig. ^a), gives rise 
to a strong coupling g. 

Fig. mb) shows the central part of current profile 
closed by a dotted ellipse in Fig. ^a.) for various w/wq 
with fixed do/d, demonstrating a larger gap for larger 
w/wq- At the boundary (x = ±w/2), the electric poten- 
tial {l/ci)dXi{x)/dx is continuous. Fig. ^c) shows the 
currents for various do/d with fixed w/wo, which shows 
the gap also grows along with do/d. These figures show 
that the current gap S grows along with both w/wo and 
do/d, and is finally saturated to \/2. 

The coupling constant in Eq. ([T|) can be rewritten as 




(16) 



where u;° = -k /\/TcL is the frequency of the 1st harmonic 
mode of the uniform resonator as can be seen in Eq. (j6]) 
and Z = \/Tfc is the impedance of the resonator<^ In 
Fig. [3] we show the coupling constant g in the plane of 
{do/d, w/wo) with Z = ^Q^^i^ From the contour plot 
of the g at the bottom of the figure the strong coupling 
regime {g ^ fiwj!) is shown to be achievable with if and 
\/d of just several multiples of wq and l/c?o- The cou- 
pling g shows a maximum, where the coupling reaches 
ultrastrong coupling regime. 

In Fig. [5] the behaviors of the current gap 5, the res- 
onator frequency LOr, and the coupling g along the di- 
agonal line w/wo = 2do/d in the plane of {w/wq, do/d) 
in Fig. [3] are shown. Note that the frequency of the 
resonator mode wj? becomes small for large c. The in- 
crease of w/wo and do/d makes the average capacitance 
of the resonator larger, and thus the resonator frequency 

smaller. As a result, the coupling g in Eq. (|16p demon- 
strates a global maximum because ujr decreases while 5 
increases. 



B. Single qubit gate 



Single qubit gate can be performed by applying an ex- 



ternal driving mode H]j 



ea^e 



^"dt+£*ae*"dt. The total 



Hamiltonian is written as Ht = H/^q 



Hd , where 77 is 



given in the rotating wave approximation (RWA) of the 
Hamiltonian in Eq. (jl3p . By using the transformation^ 
V{-f) = e''"^-''''' with 7(t) = -(e/A,.)e-*"''* and A,. = 
LJr — Wrf, we can get the transformed Hamiltonian H = 
V^HtV - iV^V given by 



H = A,.a^a H -(j 



ig{a)a^ - a(7+) + ^cr,y (17) 




djd 



FIG. 3: Coupling constant g for the single qubit case in 
the plane of {do/ d,w /wq) with the resonator impedance Z — 
50Q. g has a maximum value corresponding to a ultrastrong 
coupling. At the bottom the contour plot of g is shown. 



with Aa ^U!a- i^d- 

This Hamiltonian can be considered as if it is derived 
from the RWA of the Hamiltonian = Ara'^ a+ — 
g{a + a^)(Ty + which is identical to the Hamilto- 

nian for the circuit QED with transmon qubit in Ref. Q 
except ay instead of . Hence the analysis of qubit op- 
eration for our current-biased flux qubit can be obtained 
from Ref. by simply replacing with ay . In the dis- 
persive regime |A| ^ 5 with A = Wq — w^, for example, 
the coupling between qubit and resonator can be elim- 
inated by introducing the transformation H = W HU 

withZi = g-ii(ata_+Qa+) 



resulting in 



n 



Artt' a -\ ^'^z 



X(a'a + -)(T^ + — CTy 



(18) 



with X = 9^ / This Hamiltonian shows the resonator 
transmission shift ±;\; depending on the qubit state. 

For the present current-biased flux qubit we can get 
much stronger coupling g > A. In this case we cannot 
use above small parameter expansion any more, and thus 



5 



should calculate exact energy eigenvalues by diagonaliz- 
ing the Hamiltonian in the RWA 

which results in the eigenvalues 



i?„ = ±y(^-j +{n+l)g^. (20) 

For weak coupling regime g A with the photon number 
n = in the resonator we have £'„ ~ ±[(a;a — L0r)/2 + 
f/^/A], where the first term is the energy for the qubit 
and resonator, and the second term corresponds to the 
frequency shift of ±x- 

III. A SCALABLE DESIGN AND TWO-QUBIT 
COUPLING 

A. Scalable design 

Figure HJa) shows the schematic diagram of a scal- 
able design for the circuit QED with current-biased flux 
qubits, where dc-SQUIDs are inserted between the qubit 
and the ground plane for switching the qubit-resonator 
coupling. $xi and 3>si are the external and switching 
flux for i-th qubit, respectively. Usually superconduct- 
ing qubits are coupled with the resonator by the voltage 



-L/2 -{L-w)/2 + X - . X, + (L-w)/2 L/2 
(a) H \ 1 , ' , 1 ' F F \ H' 




FIG. 4: (a) A schematic diagram for a scalable design of 
circuit QED with current-biased three-junctions flux qubits. 
Here we show, for example, two qubits and Xi^± is given in Ap- 
pendix A. The dc-SQUIDs between qubit and ground plane 
take the role of switching the coupling between qubit and res- 
onator. The capacitance line has the width of w. (b) Current 
profiles when there are nine qubits in the circuit of (a) for 
{do/d,w/wo) = (1, 1), (5, 10), (15, 30). The current gaps ap- 
pear at qubit sites, and grow as the width w of capacitance 
line increases. 



mode or the mutual inductance. The current-biased flux 
qubit, however, is coupled with the resonator through the 
current flowing across the capacitance between the qubit 
and the resonator. Hence, if we switch off the bias cur- 
rent by piercing a half-flux quantum into the dc-SQUID 
loop, the qubit and the resonator are decoupled from each 
other. 

In Fig. IDJb) we show the current profiles in the res- 
onator for 9 qubits. In this case we use the 10th har- 
monic mode of the resonator. Similarly to the single 
qubit case of Fig. HJa) the current gap increases at the 
qubit sites along with w/wq and do/d. However, the cou- 
pling strength is different from that of the single qubit 
case. 

If there are N qubits in the circuit, we need (N-l-l)th 
harmonic mode in the resonator. For large N we use 
higher resonator mode, and thus Wr has larger value in 
Fig. El^b). The higher mode has a shorter wave length, 
and thus larger amount of bias current flows from the 
capacitance line of width w to the qubit. As a result, 
the current gap S becomes larger for many qubits [Fig. 
[5l^a)]. Since both and S have large values for many 
qubits in the circuit, the coupling g in Eq. (jl6[) increases 
along with the number of qubits N [Fig. EI^c)]. Note 
that though the high frequency of resonator field might 
be technically challenging, the resonator frequency can 




FIG. 5: (a) The current gap S for N —1, 2, 5, 9 qubits in 
the circuit of Fig. |4]as a function of w/wo along the diagonal 
line w/wq — 2do/d in the bottom of Fig. O increases 
along with w/wq and saturate to 1 finally, (b) The resonator 
frequency Wr decreases as w/wo increases. As A'^ increases, 
we need higher resonator mode and thus higher resonator fre- 
quency Wr- (c) The coupling constant g shows a maximum 
which is larger for more number of qubits. 
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be minimized, if we increase w/wq or dg/d [Fig. [5jb)]. 



B. Two-qubit coupling 

The universal gate in quantum computing requires a 
two-qubit gate in addition to the single qubit operation. 
In the scalable design of Fig. Hlja) the two-qubit Hamil- 
tonian is given by 



H. 



2q- 



fc=1.2 



k=l,2 



-ao-+fc).(21) 



In the dispersive regime |Afc| = \uJak ^ ^r \ ^ gk, we can 
obtain the coupling Hamiltonian 

^^int = ^ (^-^ + .gi52(o--lCT+2 + cr+lO--2) (22) 

by introducing a transformation similar to that in Ref. 

a 

When the qubit-resonator coupling g is strong, the 
Hamiltonian should be solved by exact diagonalization 
rather than by series expansion with a small parameter. 
The two-qubit Hamiltonian of Eq. (PT|) can be written 

-ffnubit2 and we introduce a 



as H; 



2q 



cavity ' 



transformation matrix 



qubit 1 ' 



^2 = g-«^(a*cr-l+ao-+l)-i^(aV_2+ciCT+2) ^24) 



in the same basis. Then the Hamiltonian H2q and the 
transformation matrix U2 can be represented by a block- 
diagonal matrix by slightly changing the order of basis. 

For simplicity, we consider nominally identical qubits, 
ujai = uja2, Qi = 32, and Lpi = </52, and then the lowest 
block involving the resonator photon number n = and 
1 in the Hamiltonian is written as 



2q 




(25) 



with the basis {|1 U),|0 U),|0 if)}, where | f) and 
I 1) are the qubit states, and |0) and |1) are the photon 
number states. Further the lowest block of the trans- 
formation matrix in the same basis can be evaluated by 
exactly summing a infinite series as 



U2 



I 



cos 
^sin^ 
^sin^ 



._sin^ __smv5> 



o2 V 



2 ^ 



COS^ ^ 



(26) 



Then we can easily check that if the condition tan 2(p = 
2-j2g/ /S. is satisfied, the transformed Hamiltonian H2q ~ 
Z^2-^2qW2 becomes block-diagonal further as follows: 



H- 



2q 



^/2g tan 





^tan(p Tftan^ 



(27) 



which describes the xy-type coupling between two states, 
|0 t^) and |0 4t)j with the coupling constant 



J 



9 

V2 



tan Lp 



(28) 



for n = 0. J can be explicitly evaluated with above 
condition and the interaction Hamiltonian is written by 



± 



25^ + ^ 



(cr_lcr+2 + o-+lcr-2), (29) 



where the sign is -I- for A > and — for A < 
because g > 0. For the two-qubit Hamiltonian with 
the transmon qubits H2q ~ oj^a^fl + X]fc=i 2 ^^^'^zk — 
'Ylik=i 2 9k{(i^ <^~k -|- (T_|_fca), the same coupling constant 
J can be obtained by using the transformation matrix 

jj _ g-^(a*ff-l-aff+l)-^(a+(T_2-aT+2) 

In Fig. [6lja) we show \J\/hjo^ in the plane of 
{do/d,w/wo) when there are only two qubits in the cir- 
cuit of Fig. m We set Wa/27r=2GHz for usual fiux qubits, 
a;°/27r=6GHz^ and Z = 50n,^ The strong two-qubit 
coupling J ^ fujj^ can be obtained with a few multiple 




FIG. 6: (a) Contour plot for the two-qubit xy-type interac- 
tion strength J/huJ^ in the plane of (do/d, 10/wo) when there 
are only two qubits in the circuit. The parameter values are 
Z = 50n, wJ?/27r=6GHz, and a;a/27r=2GHz. (b) The two- 
qubit coupling |J| along the dotted line in (a). For weak 
coupling I J| scales as g'^/IA], while in the limit of large g/|A| 
|J| approaches g/\/2. At the point w/wq ~ 38, the sign of 
A changes, which implies the change from ferromagnetic to 
antiferromagnetic two-qubit coupling. 
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values of wq and 1/do. The two-qubit coupling J also 
shows a maximum where J/hw^ ~ 2.4. 

Fig. [6l[b) shows the cut view of |J| along the dot- 
ted line in (a). For small g/A limit the condition 
tan2(p = 2V2g/A is written as « ^/2g/A. Then, the 
transformation matrix U2 in Eq. is reduced to U2 

in Eq. (j23p , and the coupling constant J = tan ip be- 
comes J « 5^ /A in accordance with Eq. (|22p . For large 
gr/A limit, J is given by J « ±g/^/2. 

These behaviors are demonstrated in Fig. IH^b). For 
small w/wq where the coupling g is weak, it is shown 
that J ~ g^/A. As w/wo grows, A increases because 
in Fig. EJb) decreases. At the point where the sign of A 
changes, g/A goes to infinity and J approaches ±g/\^. 
Hence, by changing parameters {w/wq, do/d) around this 
point the two-qubit ferromagnetic or antiferromagnetic 
xy-type interaction can be obtained selectively. 

IV. SUMMARY 

In summary, we propose a circuit QED architecture 
with current-biased flux qubits. The three-junctions flux 
qubit is coupled with the resonator by a bias current 
flowing through the capacitance between the qubit and 
the resonator. We introduce a capacitance line extended 
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ground plane, into the qubit loop. For much larger ca- 
pacitance almost current flows through the qubit, result- 
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Appendix A: N qubits 

Let's first consider the case that the number of qubit N 
is odd. Then the spatial part X{x) of the wavefunction 
is written by 



Ak^e-'^'^ + Bk^e-'^'' {xk,- <x < Xk,+) 



AN+ie^'-^^ + BN+^e-'^'-^- (i _ | < a- < |), 

(Al) 

where fcg is even and ko is odd among k = 
0,±l,±2,...,±(iV + 1). Here Xk^± = ± f, 

2^fco±i,± = 77^T='=T' a;_jv-i,- = — and a;jv+i,+ = 
The conditions for continuity of X{x) are given by 

Afc„-ie*^^''°-i'+ +Bfc„_ie-'^^''°-i'+ (A2) 

= Ak,e'^'""-^.+ + Bfc„e-'^^'=°-i-+ 

^fc,_ie''i^^'=- +Bfc^_ie-''t^-''- (A3) 

= Ak^e''-^'^"-- + Bfe,e-*^^'==~ 

with k ±1, ±2, ±A^, 1, and the boundary con- 
ditions at the end of resonator are 

A^N-ie-'^ + B-N-ie'^ = 0, (A4) 
An+ib'"^ + BN+ie-'"^ = 0. (A5) 

Eq. (jA2p can be written in terms of fcg and Eq. ()A3P 
in terms of fco. If we set fee = — fco + 1 and use the relation 
X-k,± = — a;fc.zp, Eqs. (|A2p and (|A3|) becomes 

A_fc^e-*¥--- +S-fc,e'^'^^^''- (A6) 

y4_fc„e-*^^''°-i'+ +B_fc„e'^^^--i'+ (A7) 
= A.k^+^e-'"^^"^-''^ + S_fc„+ie^'i^^^''-^.+ , 
respectively. 

A set of boundary conditions is obtained from sum of 
Eqs. (|A3| and (|A6l), Eqs. (jA2l) and (jA7|, and Eqs. (jA4| 
and (jASP as follows: 

{A^k^+Bk^e-''^'->"-+{Ak^+B^kJe''^^''^- (A8) 
{AN+i + B^N-i)e''^+{A^N-i+BN+i)e-''^=Q. (AlO) 



From the condition for continuity of (1 / Ci)dXi{x) / dx 
similar equations are also obtained as 

- ^ Je-^'*^^'=-+ ^ [Ak^+B^kM'^'"'^- 
c c 

= -^(A_fe^+i + i?fe^_i)e-''i^^'=- (All) 
c 

+ -(Afc_i+i3_fc,+i)e^'i"-^-, 
c 

-^(A_fc„+BfeJe-*'t^"^--i.++^(Afc„+i?_feJe''i^^''°-i'+ 
c c 

= -5(A_fe„+i + Bfe„_i)e-^¥-'=o-i.+ (A12) 

+ ^(Afe„_i+i?_,„+i)e^'?^^'=''-^.+ . 
c 

The boundary conditions of Eqs. (|A8P - (|A12p are repre- 
sented in terms of A-k + Bk- On the other hand, another 
set of boundary conditions can be obtained in terms of 
A^k — Bk as follows: 

{A-k-BkJe-''^^'"-^(Ak,-B.kJe''^^''- (A13) 
= (A_fc^+i -Sfc^_i)e-^'i^^'=— (Av-i-B-fe.+Oe*"^"^-- , 
(A„fc-i3fcJe-*'i^^^-i-+-(Afc-i3_fcJe''i^-^°--+(Al4) 

- (A_fe„+i-Sfe„_i)e-^'*^^''°-'+- (A.„-i-S-fc„+i)e'^'*^"^°-'+ , 
{AN+i-B^N-i)e''^- {A^N-i-BN+i)e~''^=0, (A15) 

-^(A_fe^-SfeJe-^'%^^'=--^(Afe-B_fcJe''*^^'=- 
c' c 

= -^(A_fc^+i-i?fc,_i)e-*'i^^'=- (A16) 
c 

-^(Afc^_i-B_fc^+i)e''i"-''-, 
c 

-:^(A_fe-i3feJe-'^-^°--+-:^(Afe-i?_feJe'^^^'=°-.+ 
c c 

= -^(A_fc„+i - Bfe„_i)e-'^-'=-i^+ (A17) 
c 

-:^(Afc„_i-i?_fe„+l)e^^-^-i.+ . 
c' 

As a result, there are two sets of boundary condi- 
tions of (i) Eqs. (|M)) - (Pl2l) and (ii) Eqs. (|JT3| - ((M7| . 
Each set can be treated as an independent eigenvalue 
problem. If the determinant of matrix corresponding to 
set (i) is non-zero while that for (ii) is zero, we have 
A-k = —Bk- Around the central qubit site the solution 
becomes ^0(2:) ^ ^0 sin ^x. On the contrary, if the de- 
terminant for (i) is zero while that for (ii) is non-zero, we 
have A^k = Bk and Xq{x) ~ cos For the case 

that the number of qubits N is even, a similar analysis 
can also be performed. 



